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TITLE OF THE INVENTION 

APPARATUS AND METHOD FOR EQUALIZING HISTOGRAM OF AN IMAGE 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application claims the benefit of Korean Application No. 2001-49810, filed 
August 18, 2001, in the Korean Intellectual Property Office, the disclosure of which is 
incorporated herein by reference. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

[0002] The present invention relates to a method of equalizing a histogram of an image, 
and more particularly, to a method of equalizing a histogram of an image using a Gaussian 
model so that the contrast can be adjusted. 

2. Description of the Related Art 

[0003] Inspecting a tone distribution of an image to be processed is essential to image 
preprocessing. A graph showing a tone distribution with respect to all pixels in the image, 
within an MxN pixel matrix using frequency by tones is referred to as an image histogram. 

[0004] When the image histogram inclines from a bright to a dark side or vice versa or is 
concentrated on a particular tone value, a corresponding image cannot be said to be good. 

[0005] When the histogram of an input image inclines to one side, tone conversion is 
used for obtaining an image with a high contrast. Essentially, the tone conversion converts 
a tone distribution of the input image according to predetermined conversion characteristics. 

[0006] Particularly, when information about an object desired to be recognized is 
unknown, histogram equalization is used as a method for automatically enhancing an image 
of the object. This method changes the tone distribution on the image so that the image 
histogram is even. In this method, probabilities of all tone values are equalized to be the 
same. 
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[0007] FIG. 1 A shows a conventional one-way histogram equalization apparatus, which 
does not consider relative brightness. FIG. 1 B shows a conventional two-way histogram 
equalization apparatus which considers relative brightness. 

[0008] The one-way histogram equalization apparatus shown in FIG. 1A includes a 
binary counter 102 for obtaining a cumulative distribution function (CDF) with respect to the 
input image through binary counting . The one-way equalization apparatus also includes a 
histogram equalizer 104 for increasing the contrast of a given image using the CDF. 

[0009] The binary counter 102 counts a number of pixels with each particular value in a 
given image, obtains a probability density function (PDF), and performs cumulative 
summation on the PDF to obtain a CDF. 

[0010] The histogram equalizer 104 applies the CDF to the given image by way of 
nonlinear mapping, thereby obtaining an image having an equalized histogram, that is, an 
image with an increased contrast . 

[001 1] However, it is difficult to control a degree of contrast enhancement in a 
conventional histogram equalization apparatus as shown in FIG. 1 A because the PDF and 
the CDF are determined based on the given image, and because the PDF and the CDF are 
used as they are for mapping performed for contrast enhancement. In addition, when the 
attribute of the given image is peculiar or is degraded due to noise, an undesirable image 
result is frequently obtained. 

[0012] In the one-way histogram equalization apparatus as shown in FIG. 1 A, relative 
brightness is not maintained. The brightness of the image obtained through histogram 
equalization has no relation with the brightness of the given image, so a bright scene is not 
discriminated from a dark scene in a video sequence, that is, in consecutive images. 

[0013] To overcome the above problems, the two-way histogram equalization apparatus 
shown in FIG. 1 B has been proposed. The two-way histogram equalization apparatus 
includes an image segmentation unit 110 to segment a given image into a subimage having 



a bright region and a subimage having a dark region . Binary counters 1 12 and 1 16 and 
histogram equalizers 1 14 and 1 18 are structurally and functionally the same as the one-way 
histogram equalization apparatus shown in FIG. 1 A . The two-way histogram equalization 
apparatus includes a mixer 120 to combine the results from the output from the histogram 
equalizers 114, 118. 

[0014] The two-way histogram equalization apparatus segments the image into a bright 
region and a dark region and separately performs histogram equalization on the segmented 
regions, so relative brightness is maintained. However, as in the one-way histogram 
equalization apparatus, when the attribute of a given image is peculiar or is degraded due to 
noise, an undesirable image result is obtained. 

SUMMARY OF THE INVENTION 

[0015] Various objects and advantages of the invention will be set forth in part in the 
description that follows and, in part, will be obvious from the description, or may be learned 
by practice of the invention. 

[0016] To solve the above-described problems, the present invention provides a 
histogram equalization apparatus, which allows a relative brightness of an image to be 
maintained and easily controls a degree of contrast enhancement by calculating a 
cumulative distribution function (CDF) using parameters estimated from a Gaussian 
distribution of the image and an error function with respect to a standard Gaussian 
distribution. 

[0017] It is a second object of the present invention to provide a histogram equalization 
method which allows a relative brightness of an image to be maintained and easily controls a 
degree of contrast enhancement by calculating a cumulative distribution function (CDF) 
using parameters estimated from a Gaussian distribution of the image and an error function 
with respect to a standard Gaussian distribution. 



[0018] To achieve the first and other objects of the invention, there is provided an apparatus 
equalizing a histogram using a CDF of an image. The apparatus includes a model 
parameter estimator estimating parameters, mean and variance, from the image using a 
Gaussian model ; an error function storing unit storing error function values based on a 
Gaussian distribution with a zero mean and a unit variance; a CDF calculator calculating a 
CDF using the parameters from the model parameter estimator and the error function values 
from the error function storing unit; and a histogram equalizer performing histogram 
equalization using the CDF calculated by the CDF calculator. 

[0019] Further, an apparatus is provided equalizing a histogram using a cumulative 
distribution function (CDF) of an image, the apparatus including: a model parameter 
estimator estimating parameters from the image using a Gaussian model; an error function 
storing unit storing error function values based on a Gaussian distribution a CDF calculator 
calculating a CDF using one of the error function values from the error function storing unit 
and the estimated parameters; and a histogram equalizer performing histogram equalization 
using the CDF. 

[0020] The apparatus further includes a model parameter modifier modifying the 
parameters, mean and variance, estimated by the model parameter estimator and outputting 
modified parameters to the CDF calculator. 

[0021] To achieve the second and other objects of the invention, there is provided a method 
to equalize a histogram using a CDF of an image. The method includes : estimating 
parameters from the image using a Gaussian model ; calculating a CDF using the 
parameters and an error function value ; and performing histogram equalization using the 
CDF. 

[0022] A method is provided to equalize a histogram using a cumulative distribution function 
(CDF) of an image, the method includes: estimating parameters including mean and 
variance, from the image using a Gaussian model; storing error function values based on a 
Gaussian distribution with a zero mean and a unit variance; calculating a CDF using the 



parameters and the error function values; and performing histogram equalization using the 
CDF calculated. 



[0023] The method further includes segmenting the image into subimages, each 
subimage including Gaussian distributions. Here, the CDF is calculated with respect to 
each of the subimages, and the histogram equalization is performed by applying the CDFs to 
the subimages. 

[0024] These together with other objects and advantages, which will be subsequently 
apparent, reside in the details of construction and operation as more fully hereinafter 
described and claimed, reference being had to the accompanying drawings forming a part 
thereof, wherein like numerals refer to like parts throughout. 



apparent by describing in detail preferred embodiments thereof with reference to the 
attached drawings in which: 



apparatus according to the present invention; 

FIGS. 3A shows a change in brightness by varying a mean from a given image; 

FIGS. 3B shows a change in contrast by varying a variance from the given image; 

FIG. 4 is a block diagram of a second embodiment of a histogram equalization 
apparatus according to the present invention; 

FIGS. 5A schematically shows a PDF of an original image during a two-way 
histogram equalization; 

FIGS. 5B schematically shows the PDF of a Gaussian model during the two-way 
histogram equalization; and 

FIGS. 5C schematically shows two CDFs after the two-way histogram equalization. 




BRIEF DESCRIPTION OF THE DRAWINGS 



[0025] The above objects and advantages of the present invention will become more 




FIG. 1 A shows a conventional one-way histogram equalization apparatus ; 
FIG. 1 B shows a conventional two-way histogram equalization apparatus ; 
FIG. 2 is a block diagram of a first embodiment of a histogram equalization 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

[0026] Hereinafter, embodiments of the present invention will be described in detail with 
reference to the attached drawings. 

[0027] FIG. 2 is a block diagram of a first embodiment of a histogram equalization 
apparatus according to the present invention. The histogram equalization apparatus shown 
in FIG. 2 includes a model parameter estimator 202, a model parameter modifier 204, an 
error function storing unit 206, a cumulative distribution function (CDF) calculator 208, and a 
histogram equalizer 210. 

[0028] The model parameter estimator 202 estimates two parameters, a mean and a 
variance, from a given image using a Gaussian model . The model parameter modifier 204 
modifies the parameters (i.e., the means and the variance) estimated by the model 
parameter estimator 202. This modification causes a shape of a cumulative distribution 
function (CDF), which is generated by the CDF calculator 208, to change, thereby enhancing 
a contrast of the image . 

[0029] The CDF calculator 208 calculates the CDF using an error function value based on a 
Gaussian distribution with a zero mean and a unit variance from the error function storing 
unit 206, and the modified parameters, i.e., the modified mean and variance from the model 
parameter modifier 204. In one embodiment, the error function storing unit 206 may be 
implemented in the form of a look-up table. The histogram equalizer 21 0 performs 
histogram equalization using the CDF calculated by the CDF calculator 208. 

[0030] The operations of the histogram equalization apparatus shown in FIG. 2 will now 
be described in more detail. Estimation of tones is essential to histogram equalization. 
One of simple methods of estimating tones is to use the Gaussian model. In the Gaussian 
model, it is satisfactory to estimate only two model parameters, i.e., a mean, ^ and a 
variance, a 2 . 
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[0031] The estimated probability density function (PDF) of a given image is given by 
Formula (1 ) 

P Fix) = _l e -^"° ! ...(1) 

[0032] Further, the CDF with respect to the Gaussian model is obtained by using 
Formula (2) 

£,(,) = T= L T le-^ 2/2 ° 2 dt = erf(x; M ,* 2 ) .-(2) 

Here, erf(x;ju,a 2 ) is an error function from the Gaussian distribution having 

(x;//,cr 2 ). 

[0033] In the Gaussian model, the CDF can be calculated referring to a table storing error 
function having a zero mean and a unit variance. Accordingly, a look-up table may 
substitute for binary counting and cumulative summation performed in a conventional 
histogram equalization apparatus. Mapping, h, used for enhancing image contrast is a 
function of the two parameters. 

[0034] The following relationship is well known: erf(x\ju,cr 2 ) = erf(ax-ju;Q,l) . In this 
case, erf \ax -//;0,1) is determined by the error function having a zero mean and a unit 
variance and az + // . The error function having a zero mean and a unit variance is stored 
in the look-up table, i.e., the error function storing unit 206, and an error function value can 
be obtained from the look-up table using cz + /i as an address. 

[0035] The CDF calculator 208 calculates an error function value using the following 
algorithm: 

for x = 0,255 
tt - ax + jU 
y = erf(ft;0$ 
CDF(x) = y; 

end. 



Here, x is a gradation value, and, for illustrative purposes, x ranges from 0 to 255. 

[0036] FIGS. 3A and 3B show effects of varying the parameters from the given image. 
FIG. 3A shows a change in brightness by varying the mean from the given image. FIG. 
3B shows a change in contrast by varying the variance from the given image. Particularly, 
varying the variance causes a slope of the error function to change. The enhanced contrast 
of an image changes depending on a change in the slope of the error function. 

[0037] FIG. 4 is a block diagram of a second embodiment of the histogram equalization 
apparatus according to the present invention. In the second embodiment, the histogram 
equalization apparatus is configured considering relative brightness. 

[0038] The histogram equalization apparatus shown in FIG. 4 includes a model 
parameter estimator 302, model parameter modifiers 304 and 404, an error function storing 
unit 306, CDF calculators 308 and 408, histogram equalizers 310 and 410, and an image 
segmentation unit 312. 

[0039] An operation of the histogram equalization apparatus shown in FIG. 4 includes 
operating one system including the model parameter modifier 304, the error function storing 
unit 306, the CDF calculator 308, and the histogram equalizer 310 and operating another 
system including the model parameter modifier 404, the error function storing unit 306, the 
CDF calculator 408, and the histogram equalizer 410. Each system processes a 
corresponding segmented subimage received from the image segmentation unit 312. The 
one system processes a bright image region and the other system processes a dark image 
region . The operation of each system is the same as described in FIG. 2, with the 
exception that each system processes different image regions. 

[0040] The operations of the histogram equalization apparatus shown in FIG. 4 will now be 
described in detail. FIGS. 5A through 5C schematically show estimation of tones during a 
two-way histogram equalization process. FIG. 5A shows a PDF of an original image. 
FIG. 5B shows the PDF of the Gaussian model. FIG. 5C shows two CDFs after the 
two-way histogram equalization process is performed. 
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[0041] Estimated model parameters are applied to two Gaussian distributions separately. It 
is assumed that two vectors are F = (F 09 F ] ,.„F MN _ J ) t and C = (C 09 C 19 ...,C m _ ] ) t . The 
real parts of F and Care expressed by/ = (f 0 ,f I ,-,fMN-i) t and c-ic^Cj^c^) 1 , 
respectively, and can be expressed by an MxN-matrix having f{m,n) and c(m,n) as elements. 

[0042] Cluster information c(rr?,n) indicates a distribution to which a pixel f{m,n) belongs. 
Each pixel Km,n) has a value of x e [0,L] . The PDF of a given image may be expressed 
by Formula (3). 

P F =t 4 P[f = x\c = mc = i] -(3) 

Here, P[f = x\c = i] denotes a Gaussian distribution, and P[c = i] denotes a class 
Q probability. 

ffl [0043] A given image f is considered as incomplete data losing cluster information c. 



hi 



Estimation of tones with respect to complete data (f, c) is expressed by Formula (4). 

0 = argmaxP[/,c|0] ...(4) 
Here, 0 = (p 29 aj ,fi 2 ,a62 2 ,pc } y . 



W [0044] According to AP. Dempster, N.M Laird, and D.B. Rubin, "Maximum Likelihood 

m 

Q from Incomplete Data via the EM algorithm," J. Royal Statistical Soc. Ser. B,Vol. 39, 

* y n9.1 , pp. 138, 1977 ] and R.J.A. Little and D.R. Rubin, Statistical Analysis with Missing 

Data, John Wiley & Sons, New York, 1987 ], an expectation maximization (EM) algorithm 

provides a maximum likelihood for Formula (4). 

[0045] The EM algorithm is accomplished by repeating the following two steps : 
E-step;Q(0\d T ) =E[logP[f,c\0][f,6*] ...(5) 

M-step;0 T+7 -arg maX e(00 r ) ...(6) 
0 
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[0046] Cluster information is determined in the M-step, and parameters for each cluster 
are individually obtained. The cluster information, c, can be obtained through hypothesis 
testing for individual elements. 



c = < 



J - if Plf\c = IF ]p T d > p[f\c = 2 9 0 r ](l - p T cl) 
2. otherwise 



■(7) 



[0047] After each of all pixels is classed as one of two clusters according to the 
hypothesis testing, the mean and the variance are individually calculated with respect to 
each cluster. 

[0048] After all parameters are estimated according to the EM algorithm, the CDFs of two 
Gaussian distributions are determined using the following two formulas : 

C F ,(x) = erf(x;M 1 ,<7j) -( 8 ) 

C F2 (x) = erf(x;ju 2 ,a 2 ) ...(9) 

The CDFs are used as mappings h 1 and h 2 in a two-way histogram equalization. 
[0049] The given image is segmented into two subimages f and f by Formula (1 0). 



ti, x \f(m,n),if c{m,n) = i 

f {m,n) = \ . ...(10) 

[0, otherwise 



C F1 and C F2 are applied to the two subimages f* and ^as the mappings 
h' :[0,L c ]^[0,L c ] and h 2 :[L C> L]-+[L C ,L] , respectively. 



h'w')=- —~ c F1 (f) ...(11) 

C F1 (L c )-C F1 (0) 
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h\f 2 ) = -~ ^ C F2 {f) -(12) 

C F2 (L c )-C F2 (0) 



Here, L c is a value satisfying Formula (13). 

P[f = L c \c = l]pc 1 =F[f = L c \c = 2\l-pc 2 ) ...(13) 

[0050] If it is assumed that the enhanced subimages have a uniform distribution, an 

. , x* u u ■ ■ (0+L c )pcl (L c +L)(l-pc2) 
average bnghtness of the enhanced entire image is + — . Even 

2 2 

after applying the two-way histogram equalization, the average brightness of the enhanced 
|5j entire image is still a function of L c , and changes in the relative brightness between scenes 



•3P are maintained. 



[0051] In the histogram equalization apparatus according to the present invention, a 
binary counter and a cumulative summation unit required by a conventional histogram 
equalization unit can be replaced by the error function storing unit realized as a look-up 



W table, thereby simplifying the structure and allowing the result of histogram equalization to be 



easily controlled using model parameters. 



■pes 

% 

flJ [0052] In addition, because the histogram equalization apparatus according to the present 

invention controls the parameters used to obtain the CDF using the model parameter 
modifier, it is easy to control the degree of contrast enhancement. 

[0053] The histogram equalization apparatus according to the present invention can 
perform histogram equalization without noise or pattern amplification through model 
parameter modification. Moreover, when the histogram equalization apparatus according to 
the present invention is applied to a two-way histogram equalization, relative brightness of 
given images is maintained, thereby the present invention can be easily applied to a video 
sequence. 
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[0054] While this invention has been particularly shown and described with reference to 
preferred embodiments thereof, it will be understood by those skilled in the art that various 
changes in form and details may be made therein without departing from the spirit and scope 
of the invention as defined by the appended claims. 
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